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Abstract. This paper is based on a formulation of the Navier-Stokes equa- 
tions developed by P. Constantin and the first author (arxiv:math. PR/0511067, 
to appear), where the velocity field of a viscous incompressible fluid is written 
as the expected value of a stochastic process. In this paper, we take N copies 
of the above process (each based on independent Wiener processes), and re- 
place the expected value with times the sum over these N copies. (We 
remark that our formulation requires one to keep track of N stochastic flows 
of diffeomorphisms, and not just the motion of N particles.) 

We prove that in two dimensions, this system of interacting diffeomor- 
phisms has (time) global solutions with initial data in the space C 1,Q which 
consists of differentiablc functions whose first derivative is a Holder continuous 
(see Section [3] for the precise definition). Further, wc show that as N — > oo 
the system converges to the solution of Navier-Stokes equations on any finite 
interval [0, T] . However for fixed TV, we prove that this system retains roughly 
O(jf) times its original energy as t — > oo. Hence the limit N — > oo and T — > oo 
do not commute. For general flows, we only provide a lower bound to this ef- 
fect. In the special case of shear flows, we compute the behaviour as t — » oo 
explicitly. 



1. Introduction 

The Navier-Stokes equations 

(1.1) d t u + (u ■ V)ti - vl\u + Vp = 

(1.2) V-u = 

describe the evolution of a velocity field of an incompressible fluid with kinematic 
viscosity v > 0. These equations have been used to model numerous physical 
problems, for example air flow around an airplane wing, ocean currents and meteo- 
rological phenomena to name a few p||15U18) . The mathematical theory (existence 
and regularity [3[T4]) of these equations have been extensively studied and is still 
one of the outstanding open problems in modern PDE's [HIE]. 

The questions addressed in this paper are motivated by a formalism of (jl.ljl 
(|1.2p developed in 9: (equations (|2.1|) - (12.2p below). This formalism essentially 
superimposes Brownian motion onto particle trajectories, and then averages with 
respect to the Wiener measure. In this paper, we take N independent copies of 
the Wiener process and replace the expected value in the above formalism with 
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times the sum over these N independent copies (see equations (|2.5p - (|2.6[) for 
the exact details). In the original formulation the random trajectory of a particle 
induced by a single Brownian motion interacts with its own law. This is essentially 
a self-consistent, mean-field interaction. In this paper, we replace this with N copies 
or replica whose average is used to approximate the interaction with the processes 
own law. This technique has been extensively used in numerical computation (e.g. 
[21[T7l[T9] ) . We remark that in our formulation, we are required to keep track of N 
stochastic flows of diffeomorphisms, and not just the motion of N different particles, 
as is the conventional approach. 

We study both the behaviour as N — + oo and t — > oo of the system obtained. The 
behaviour as N — > oo is as expected: In two dimensions on any finite time interval 
[0,T], the system converges as N — > oo to the solution of the true Navier-Stokes 
equations at rate roughly 0(-J=). In three dimensions, we can only guarantee this if 
we have certain apriori bounds on the solution fTheorem l4.1[) . These apriori bounds 
are of course guaranteed for short time, but are unknown (in the 3-dimensional 
setting) for long time [3[8] . 



(a) N = 2 
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(b) N = 8 



Figure 1. Graph of ||wf]|f,2 vs time 



for fixed N is less intuitive. For the 



2 L 2 vs time, with 



At first glance, the behaviour as t — 
2-dimensional problem, Figures 1(a) and 1(b) show a graph of 

N = 2 and N — 8 respectively. A little reflection shows that this behavior is not 
completely surprising. The dissipation occurs through the averaging of different 
copies of the flow. With only N copies, one can only produce dissipation of order 
1/N of the original energy It is tempting, to speak of the inability to represent 
the correct interaction of small scale structures with such small number of data. 
However, we can not make this precise and since each of the objects being averaged 
is an entire diffeomorphism with an infinite amount of information it is unclear 
what this means. 



These computations were done using a 24 X 24 mesh on the periodic box with side length 
2ir. The initial vorticity was randomly chosen, and normalized with ||uo||r2 = 1. The behaviour 
depicted in these two figures is however characteristic, and insensitive to changes of the mesh size, 
length, or diffusion coefficient 



A STOCHASTIC-LAGRANGIAN PARTICLE SYSTEM FOR NAVIER-STOKES 



3 



In Section [5] we obtain a sharp lower bound to this effect. We show (Theorem 
EUD that 

limsupE||Vut||^2 ^ \\u \\ 2 L 2 , 

where L is a length scale. Further, we explicitly compute the t — > oo behaviour in 
the special case of shear flows and verify that our lower bound is sharp. 

We remark that we considered the analogue of the system above for the one 
dimensional Burgers equations. As is well known the viscous Burgers equations have 
global strong solutions. However preliminary numerical simulations show that the 
system above forms shocks almost surely even for very large N. We are currently 
working on understanding how to continue this system past these shocks, in a 
manner analogous to the entropy solutions for the inviscid Burgers equations, and 
studying its behaviour as t — > oo and N — ► oo. 

We do not propose this particle system as an efficient particle method for nu- 
merical computation. Though there may be special cases were it may be useful, in 
general the computational cost of representing N entire diffcomorphisms is large. 
Rather we sec it as an interesting and novel rcgularization which might give useful 
insight into the structure and role of dissipation in the system. 

2. The particle system 

In this section we construct a particle system for the Navier-Stokes equations 
based on stochastic Lagrangian trajectories. We begin by describing a stochastic 
Lagrangian formulation of the Navier-Stokes equations developed in [SHU]. 

Let W be a standard 2 or 3-dimensional Brownian motion, and uq some given 
divergence free C 2 ' a initial data. Let E denote the expected value with respect 
to the Wiener measure and P be the Leray-Hodge projection onto divergence free 
vector fields. Consider the system of equations 

(2.1) dX t (x) = ut(X t (x)) dt + V2v dW t , X (x) = x, 

(2.2) ut = EP[(V*n)(«ooy t )]. Y t =Xt 1 . 

With a slight abuse of notation, we denote by X t the map from initial conditions 
to the value at time t. Hence X t is a stochastic flow of diffcomorphisms with Xq 
equal to the identity and Y t the "spatial" inverse. In other words, Y t : X t (x) i— > x. 
Also by V*Yt we mean the transpose of the Jacobian of map Y t . Observe that 
(V*Y t )(uo o Yt) can be viewed as a function of x where both the Jacobian and 
the vector field uq o Y t to which it is applied are both evaluated at x. We impose 
periodic boundary conditions on the displacement Xt(y) = X t (y) — y, and on the 
Leray-Hodge projection P. 

In [9l[TT] it was shown that the system (|2.ip - (|2.2p is equivalent to the Navier- 
Stokes equations in the following sense: If the initial data is regular (C 2 ' Q ), then the 
pair X : it is a solution to the system (|2.ip - (|2.2p if and only if u is a (classical) solution 
to the incompressible Navier-Stokes equations with periodic boundary conditions 
and initial data uq. 

We digress briefly and comment on the physical significance of (|2.ip - (|2.2p . Note 
first that equation (|2.2|) is algebraically equivalent to the equations 

(2.3) u t = -A _1 V x uj t 

(2.4) u t =El(VX t )u ]oY t 
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This follows by direct computation, and was shown [6] and [9] for instance. We 
recall that (|2.4|) is the usual vorticity transport equation for the Euler equations, 
and (|2.3p is just the Biot-Savart formula. 

Thus in particular when v = 0, the system (|2. 1|) — (|2.2[) is exactly the incompress- 
ible Euler equations. Hence the system (|2.1|) ~ (|2.2p essentially does the following: 
We add Brownian motion to Lagrangian trajectories. Then recover the velocity u 
in the same manner as for the Euler equations, but additionally average out the 
noise. 

We remark that the system (j2.ip ~ (|2.2p is non-linear in the sense of McKean 
[20] . The drift of the flow X depends on its distribution. However in this case, 
the law of X alone is not enough to compute the drift u. This is because of the 
presence of the V*Y term in (12. 2p . which requires knowledge of spatial covariances, 
in addition to the law of X. In other words, one needs that law of the entire flow 
of diffeomorphism and not just the law of the one-point motions. 

We now motivate our particle system. For the formulation (I2.ip - (|2.2p above, 
the natural numerical scheme would be to use the law of large numbers to compute 
the expected value. Let (W l ) be a sequence of independent Wiener processes, and 
consider the system 

(2.5) dX}> N = u? (XI> N ) dt + \f2v dWl Yt N = (x^ N ) ' 

1 N 

(2.6) u? = ^p[(v*y;> n )(u oy;< n ) 



with initial data Xo(x) — x. We impose again periodic boundary conditions on the 
initial data uo, the displacement \t(x) = X t {x) — x, and the Leray-Hodge projection 
P. 

We remark that the algebraic equivalence of (12.21) and (|2.3p - (l2.4p is still valid 
in this setting. Thus the system (j2 . 5|> — (j2 . 6(1 could equivalently be formulated by 
replacing equation (|2.6p with the more familiar equations 

1 N 

(2-7) w f = i^[(V^> 

i=l 



and 

(2.8) = -A" X V x w f . 

Finally we clarify our previous remark, stating that the above formulation re- 
quires us to keep track of N stochastic flows, and the knowledge of the one point 
motions of X\' N alone is not sufficient. The standard method of obtaining a solu- 
tion to the heat equation (assuming the drift u is time independent) would be to 
consider the process (|2.ip . and read off the solution 9 by 

e t {a)=W a oX t {a) 

where 8q is the given initial temperature distribution. Thus knowing the trajecto- 
ries (and distribution) of the process X starting at one particular point a will be 
sufficient to determine the solution 9 t at that point. 

This however is not the case for our representation. The reason is twofold: 
First, the representation (|2.ip - (|2.2p involves a non-local singular integral operator. 
Second our representation involves composing with the spatial inverse of the flow 
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Xt, and then averaging. If we for a moment ignore the non-locality of the Leray- 
Hodge projection, determining ut at one fixed point a one would need the law of 
Yt (a) , for which the knowledge of X t (a) alone is not enough. One needs the entire 
(spatial) map X t to compute the spatial inverse Y t (a). 

The above is not a serious impediment to a numerical implementation. Given an 
initial mesh A, we first compute X\ ,N on this mesh. By definition of the inverse, one 
knows Y%' on the (non-uniform) mesh X l t ' (A), after which one can interpolate 
and find Y^' N on the mesh A. In two spatial dimensions, global existence and 
regularity (Theorem 13. 5p together with incompressibility will show that this mesh 
does not degenerate in finite time. 

This surprisingly is not the case for the (one dimensional) Burgers equations. 
Numerical computations indicate the mesh A almost surely degenerates in finite 
time for non-monotone initial data, and the solution 'shocks' almost surely. Thus, 
while the system (|2.5I) ~ (|2.6I) appears natural, and convergence as N — ► oo is to 
be expected, caution is to be exercised. We suspect that the results (existence, 
convergence, etc.) proven in this paper for the system (|2.5|) - ([2.6|) are in fact false 
for the Burgers equations. This is indeed puzzling as global existence, and regularity 
for the viscous Burgers equations is well known. It further underlines the fact that 
the finite N approximation modifies the dissipation in a different way then other 
approximation such as a spectral approximation. 

In the next section, we show the existence of global solutions to (|2.5p ~ (|2.6p in 
two dimensions. In section^ we show that the solution to (|2.5p - (|2.6[) converges 
to the solution of the Navier-Stokes equations as N — + oo. Finally in Section [5] we 
study the behaviour of the system ()2.5p - (l2.6p as t — > oo (for fixed N), and partially 



explain the behaviour shown in Figures 1(a) and 1(b) 



3. Global existence of the particle system in two dimensions. 

In this section we prove that the particle system (|2.5p ~ (|2.6p has global solutions 
in two dimensions. Once we are guaranteed global in time solutions, we are able 
to study the behaviour as t — > oo, which we do in Section [5] We remark also that 
as a consequence of Theorem 13.51 (proved here), our convergence result as N — > oo 
(Theorem 14. ip applies on any finite time interval [0, T] in the two dimensional 
situation. 

We first establish some notational convention: We let L > be a length scale, 
and assume work with the spatial domain is [0,£] d , where d ^ 2 is the spatial 
dimension. We define the non-dimensional LP and Holder norms by 




NL,a= E Lk \ Dm <+ E L lml \\D m M 

\m\=k |m|<fc 
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We denote the Holder space C k,a and Lebesgue space LP to be the space of functions 
u which are periodic on [0, L] d with < oo and < oo respectively. 

Let W 1 , . . . , W N be N independent (2 dimensional) Wiener processes, with fil- 
tration Tt- 

Proposition 3.1. Let to 0, and u} () , . . . , uf g be T ta -measurable, periodic mean 
functions such that the norms \\u\ o \\i^ a are almost surely bounded. Then there exist 
T = T(a, ||wt ||i,a, ■ ■ ■ , || w tllli,a) such that the system 



(3.1) X\ Q , t {x) = x + u s (Xt 0tS (x)) ds + V^(W? - Wi ), Yi 0tt = (XIJ- 1 

J t 

1 N 

(3.2) Ut (x) — [(v*i? 0i *) « o Yi 0tt )] (x) 



N 

1=1 

has a solution on the interval [to , T] . Further there exists deterministic U = 
U(a, \\ul \\i,a) such that 

( 3 - 3 ) SU P ll p [(V**? ,tX o1 ?o 1 t]|| lo <^ 

te[t ,T] 

almost surely. Consequently, u £ C(\to, T], C '"), and ||u||i )Ct ^ U. 

Proposition 13. II is proved in Appendix [Al 

Definition 3.2. We call X 1 , u in Proposition ^. II the solution of the system (|3.ip — 
— (|3.2p with initial data u\ Q , . . . , . 

The functions X 1 ,Y % are not periodic themselves, but have periodic displace- 
ments: Namely, if we define 

(3-4) Xi(y)=Xi(y)-y 

(3.5) 4(x)=Yi(x)-x, 

then X 1 are periodic. 

We remark that if to = and the ^J o 's are all equal, then the system ()3.1[) — 
(|3.2[) reduces to (|2.5p ~ (|2.6p . However when formulated as above, solutions can be 
continued past time T by restarting the flows X % , as in the Lemma below. 

Lemma 3.3. Say to < t\ < t-z, and X\ s> u s solve (|3.1[) - (|3.2[l on [£0^2], with 
initial data u\ Q . For any t > to define 

t4 = p[(v*it*K° 

LetX^ „. solve (|3.ip ~ (l3.2p on [ti^tz] with initial data u\ t . Then for all s S [£1,^2] 
we have u l s = u\ and 

X t ,s = X t t ,s ° X t M 

almost surely. 

Proof. The proof is identical to the proof of Proposition 3.3.1 in [TT], and we do 
not provide it here. □ 

For the remainder of this section, we assume without loss of generality that 
to = (we allow of course To to be non-trivial). For ease of notation, we use X s 
to denote Xo }3 . We now prove that the system (|3.ip - (|3.2p has global solutions in 
two dimensions. This essentially follows from a Beale-Kato-Majda type condition 
lj, and the two dimensional vorticity transport. 
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Lemma 3.4. If u is divergence free and periodic in M. d , then for any a £ (0,1), 
there exists a constant c — c(a, d) such that 



|Vu||-. <c|M|--, ( 1 + ln 4 



where uo = V x u. 



The lemma is a standard fact about singular integral operators, and we provide 
a proof in Appendix [Bl for completeness. 



Theorem 3.5. In 2 dimensions, the system (j3.ip ~ p.2p has time global solutions 
provided the initial data periodic, T§ -measurable with ||w lli,a bounded 

almost surely. In particular we have global solutions to (|3.ip — 13 . 2(1 in two dimen- 
sions, if Uq = ■ ■ ■ = Uq = wo is deterministic, Holder 1 + a and periodic. 

Proof. Taking the curl of (|3.2p gives the familiar Cauchy formula [61I9HT0] 

1 N 

(3.6) c^ = -]T[(v^K]or t \ 

»=i 

where uj t = V x w t . In two dimensions reduces to 

1 N 

(3.7) wt = -Y,"o<>Yt- 

i=l 

Taking Holder norms gives 

JV 

(3-8) INIa<^EH W -IL( 1 + ll VY< ll"-) ^ 

1=1 

Now differentiating (|3.ip gives 

VX l t =I + I {Vu s o XDVXlds a.s. 
Jo 

Taking the L°° norm, and applying Gronwall's Lemma shows 

||VA^|| iOC < exp J ||Vu s ||--! ds^j a.s. 

Recall V • u = 0, and hence det(VAT l ) = 1 almost surely. Thus the entries of VY 
are a polynomial (of degree 1) in the entries of VX. This immediately gives 

(3.9) \\VYt\\ Lm ,^exp(cJ ||Vu-||-.<- d* 

almost surely. Combining this with (|3.8p gives us the apriori bound 

N 



( C 

y_ ii w oL ex p c / ii 

i=l V ^0 



Applying Lemma 13.41 gives us 
||Vu|| L oc < c||o;|| L oc ( 1 + ln 



Ml; 



L°° 

K 

1 11/ 1» I 

iV Z^i=l H'-oIlL- 
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Note that the function x ln + - is bounded, so the last term on the right can be 
bounded above by some constant cq- For the remainder of the proof, we let cq = 
Co ( 1 1 ojq 1 1 q j o) denote a constant (with dimensions that of u) which changes from line 
to line. Thus 




and hence 



Vu t \\ L oc ^1 + J || Vitals ds\ < c + c||w t || L c 



Integrating gives us the apriori bound 

J || Vu t \\ L oc < exp (c t + cj \\u) t \\ L <*>\ - 1 

(3.11) s$ c Q te Cot 

since (JXTJ) implies \\uj t \\ L <~ < j^O^olU 00 - 

Now if we set io\ — luq o Y t l , then f|3. 9|) and the apriori bound (|3.1ip gives 

(3.12) \\col\\ a H^l^exp^te 00 *) 
If u\ is as in Lemma 13.31 then ()3.12p shows 

(3.13) ||Vt4|| Q ^c\\4\\ a <c 1^0 IU exp (c i e co *). 

Since the mean velocity is a conserved quantity, a bound on ||VwJ|| Q immediately 
gives a bound on ||ttj||i a , which in conjunction with local existence (Proposition 
13. ip . and Lemma [3 . 3 1 concludes the proof. □ 

4. Convergence as N — > oo 

In this section, we fix a time interval [0, T], and show that the particle system 
(|2.5p ~ (|2.6p converges to the solution to the Navier-Stokes equations as N — > oo. 
The rate of convergence is 0(-i=), which is comparable to the convergence rate of 
the random vortex method [3lll6j. As mentioned earlier, the system is intrinsically 
non-local, and propagation of chaos [20] type estimates are not easy to obtain. 
Consequently convergence results based on spatially averaged norms are easier to 
obtain, and we present one such result in this section, under assumptions which are 
immediately guaranteed by local existence. 

Theorem 4.1. For each i,N, let X l,N , u N be a solution to the particle system 
(|2.5p - (|2.6|) with initial data uq on some time interval [0, T]. Let u be a solution 
to the Navier-Stokes equations (with the same initial data) on the interval [0,T]. 
Suppose U is such that 

sup ||Vut|| L 2 ^ — and sup 
te[o,T] L te[o,T] 

Then (u N ) — > u in the following sense: 



T7 



u 

< — 
L 



lim sup E \\uf — Ut || r 2 = 
N ->°°te[o,T] 
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We remark that given C 1,a initial data, local existence (Proposition l3.1j) guaran- 
tees that the conditions of this theorem are satisfied on some small interval [0, T]. In 
two dimensions, Theorem 13 . 5 1 shows that the conditions of this theorem are satisfied 
on any interval finite [0, T]. 

The proof will follow almost immediately from the following Lemma. 



Lemma 4.2. Let v% N = P[(V*Y t 1 '" )u o Y t 1 '"} be the 



i,N\ 



ri.N-i 



:th 



summana in 



2~6l). Then 



satisfies the SPDE 



(4.1) d U \> N + (nf-VK 



uAu t ' + (V u t )u t 



N 



T7 i,N 



dt- 



2vVui' N dWi = 



PT N 

z=l 



and u N satisfies the SPDE 

(4.2) duf + [(uf • V)uf - i/Attf + Vpf ] dt + 



Remark. We draw attention to the fact that the pressure term in (|4. 2|) has bounded 
variation in time. 

Proof. We first recall a fact from [HUH] (see also [H]). If X is the stochastic flow 

dX t =u t dt+ V2v dW t 

and Yt — X^ 1 is the spatial inverse. Then the process 0t = f(Y t ) satisfies the 
SPDE 

(4.3) 



d6 t + (lit • V)0 4 df - vt\B t dt + V2vV6 t dW t = 



This immediately shows that Y^' N and v\' N 



u o Y t l ' N both satisfy the SPDE 
()4.3p . For notational convenience, we momentarily drop the iV as a superscript and 
use the notation v 1 ' 3 to denote the j th component of v % . 
Now we set w\ = (V*y/)wJ and apply Ito's formula: 



dw t ' 3 



d{d 3 Y?) ■ v\ + (a,i?) • d^ + d (a,-*?-*, ^< fc ) 

(fyl?) • [-(rit • V)uj + ^Awj] dt - V^djYj ■ (Vvt dWl) + 
+ v\ ■ [-{(djUt) ■ V)Y t l - (u t ■ V)d 3 Y; + vAdjYj] dt- 
- V^v] ■ (VdjYj dW l t ) + 2vd] l Y t i ' k d l v i t ' k dt 

\-{u t ■ V)w\ + vAw\ - {V*u t ) ■ wl] dt - V2^Vw\ dW\ 



Restoring the dependence on N to our notation, since u 



i,N 



Pw I,,v we know that 



i.N 

H 



i.N 



i,N 



for some function q l t ' N . Thus 



ofit t = dw t 



dt 



2vVw\' N dW l t 



dt - V2^Vui' N dW? 



+ 



( U f ■ vxvgr ) + ) - (v*o(v«r v ) 



i.iV 



i,N \ 



dt- 
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'2i/V(V#")dW? +d(V?i- 



If we define P 4 4,JV by 



Jo Jo 



5/ 



then we have 

(4.4) dtff + Uu» ■ V)u\' N - uAul' N + (y*u?W 



dt 



+ d(S/Pt' N ) + \/2vVu l ^dWl = 0. 

Notice that u l,N is divergence free by definition, and thus Vu ,,JV dW% is also diver- 
gence free. Thus d(V P*' ), the only other term with possibly non-zero quadratic 
variation, must have a divergence free martingale part. Since the martingale part 
of d(VP l ' N ) is also a gradient, it must be 0. Thus 

d(vpt N )=v P i' N dt 

for some function p\' N , which proves (|4.1[) . 

The identity (|4. 2[) now follows by summing (|4.1[) in i, dividing by N, and defining 

Pt by 

AT 



1. 



AM 



1 E 

7Y f-* 



i.iV 
Pi 



□ 



Proof of Theorem \4-l\ Let u be a solution of the Navier-Stokes equations, with 



initial data Mo> and set v N = u N 



Then / satisfies the SPDE 



dvf + Of • V)u t eft + (ut • V)«f (it - i/Auf + V(pf - p t ) dt+ 

N 



>2v 



N 



i=l 



where p is the pressure in the Navier-Stokes equations, and p the pressure term 
in (Oil . 

Thus by Ito's formula, 



id\ 



D N \\ 



+ «,(< • V)u t ) dt + v\\\/v 



th? dt + 



2v 
~N 



N 



N 



f/f. 



i=l i=l 

Here the notation (/, g) denotes the L 2 innerproduct of / and g. Taking expected 
values gives us 

and by Gronwall's lemma we have 



E U - N 



2 2vU 2 ut 

V't' II r2 ^ -TTTTTte L 



concluding the proof. 



□ 
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5. Convergence as t — ► oo 

In this section, we fix N, and consider the behaviour of the system (|2.5[) — (|2.6() 
as t — » oo. We show that the system (|2.5p - (|2.6p does not dissipate all its energy as 
t — ► oo. Roughly speaking we show 

UmsupE||Vu t ||i2 >0(±), 

t — >oo 

with constants independent of viscosity. This is in contrast to the true (unforced) 
Navier-Stokes equations, which dissipate all of its energy as t — > oo (provided of 
course the solutions are defined globally in time). 

In general we are unable to compute exact asymptotic behaviour of the system 
(|2.5|) - (|2.6|) as t — > oo. But in the special case of shear flows, we compute this 
exactly, and show that the system eventually converges to a constant, retaining 
exactly times its initial energy. 

For the remainder of this section, we pick a fixed N £ N and for notational conve- 
nience we omit the superscript N. We begin by computing exactly the asymptotic 
behaviour of the system (|2.5p ~ (|2.6p in the special case of shear flows. 

Proposition 5.1. Suppose the initial data Uq(x) = (())o(x 2 ),0) for some C ' a pe- 
riodic function <j>. If u is the velocity field that solves the system (|2.5p - (|2.6p with 
initial data Uo, then 

(5.1) km Euj t (x) 2 = 1 Hwolfe 

where lu = V x u is the vorticity. 

Proof. Let X 1 , Y l be the flows in the system (|2.5p ~ (|2.6p . and as before define u l to 
be the i th summand in (|2.6|) . and oj' l = ojqoY' 1 . 

First note that the SPDE's for u l and u (equations (14. ip and (|4.2p ) are all 
translation invariant. Thus since the initial data is independent of x\, the same 
must be true for all time. Since it 1 , u are divergence free, the second coordinate 
must necessarily be 0, and the form of the initial data is preserved. Namely, 

u t (x) = (0 t (x 2 ),O) and ui(x) = (4>l(x 2 ), 0) 

for some C([0, oo), C 1,Q ) periodic functions (j> l ,<j>. 
Now the definition of X 1 shows that 

\y 2 + V2vW t J 

and hence 

Y i ( ^ _ ( x l +^ 1 {x 2 ) \ 

Yt{x) -\x 2 -v^w:' 2 )- 

Recall A 1 , \i % are as in (|3.4p . (|3.5p . and here the notation A 2,1 to denotes the first 
coordinate of A\ This immediately shows 

lo\{x) =u)ooY t l (x) 

= -d 2 M*2 - V^wi* 2 ) 

where W 1,2 again denotes the second coordinate of the Brownian motion W . 
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Now using standard mixing properties of Brownian motion [131 Section 1.3], (or 
explicitly computing in this case) we know that for every x £ [0,£] 2 

2 



(5.2) 
and 



lim EuJ l t (xY 



lim E 



-d 2 (f> (x 2 



\\vofi? 



lim E u>l (x)ujI (x) 

t — >oo 



lim E 

t— >oo 



— I d 2 

[0,i] = 

(5.3) 
when i ^ j. 

Now by two dimensional Cauchy formula ([3.7 




d 2 4>o(x 2 
E<9 2 <Ma>2 - 



2vwl' 2 ) m 2 M*2 - 



'2vW 3 t ' 2 ) 



i 



since m our case, 



U3i 



,N 



N 



N ^ 

»=i 



uj o Y< 



loq). Thus 



w JV i-l 

(5.4) E^(x) 2 = -^E W „o Y t \x) 2 + T« E E ° ° Y * ( x 



N 2 

i=l i=2 j=l 

the proof is complete. 



□ 



N 2 

and using ([5.2p and ([5 

We remark that all we need for (j5.1[) to hold is the identities (j5.2[) and (|5.3p . 
Equality (|5.2[) is guaranteed provided reasonable ergodic properties of the flow X\ 
are known. Equality (|5.3() is guaranteed provided the flows X\ and Jff eventually 
decorrelate. 

While we are unable to guarantee these properties for a more general class of 
flows, we conclude this section by proving a weaker version of (jS.ip for two dimen- 
sional flows with general initial data. 

Theorem 5.2. Let X l ,u be a solution to the system (j2.5|) - ([2.6p with (spatial) 
mean zero initial data uo € C 1,Q and periodic boundary conditions. Suppose further 
u G C([0,oo),C 1 ' '). Then 



(5.5) 



lim sup E || Vut|| r2 ^ 



1 



ii 2 



t^oo NL 2 

Note that the assumption u € C([0, oo), C 1,Q ) is satisfied in the two dimensional 
situation with C 1,a initial data (Theorem 13. 5p . The proof we provide below will 
also work in the three dimensional situation, as long as global existence and well- 
posedness of (|2.5p - (|2.6p is known. 

As is standard with the Navier-Stokes equations, the condition that uq is (spa- 
tially) mean zero is not a restriction. By changing coordinates to a frame moving 
with the mean of the initial velocity, we can arrange that the initial data (in the 
new frame) has spatial mean 0. 

Finally we remark that the lower bound in inequality (|5.5p is sharp, since in the 
special case of shear flows we have the equality (15. ip . However we are unable to 
obtain a bound on liminf E|| Vu t || 2 2 . 
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Proof of Theorem [5H As before let u\ = P [(V*Yj)uo ° Yj] . Using Lemma l4~2l 
and Ito's formula we have 



(5.6) id | 



Note that 



[(u • V)u l - v/\u l + (V*uK + Vp l ] d£+ 



+ V2i/u 4 • (Vm* dIU 1 ) = 1/ Vm 1 dt 



u* ■ ((V*uK) = / ((VuX) • u 



u l • (u ■ V)u l = 



Thus integrating (|5.6j) in space, and using the fact that u l is divergence free gives 

,2 



and hence 



H\\L 2 



|mq|| L 2 almost surely. 



Now suppose that for some e > 0, there exists to such that for all t > to 

'|Vut||/2 < ttW \\uq\\ 2 L 2 - e. 



NL 2 



Using Ito's formula and (|4.2| gives 



JV 



~d M 2 + u-[(u- V)u - uAu + Vp] dt + u ■ (Vh 1 dW l ) 



N 



d,i 



Integrating in space, and taking expected values gives 

N 



-^d t E||u t || 2 2 =E 



ui\\ L 2 ~ ||Vu t || 2 2 



1 N 



N 2 L 2 ^— ' " 

i=l 



iVutl 



1 



1 * 
272 Z)H Uo 



| 2 2 -E||V Ut || 2 2 



ll u o||?,2 - E llVutllxa 



for t ^ to. Here we used the Poincare inequality to obtain the second inequality 
above. Note that we have assumed that the initial data has (spatial) mean 0. Since 
the (spatial) mean is conserved by the system (|2.5l) - (|2.6p , Ut also has (spatial) mean 
zero, and our application of the Poincare inequality is valid. 

Now, the above inequality immediately implies E||ii t || 2 2 becomes arbitrarily large 
as t — * oo. This is a contradiction because 

1 N 1 N 



\\ u t\\ L 2 
holds almost surely. 



i=l 



a 
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Appendix A. Local existence. 



In this appendix we provide the proof of Proposition 13.11 A similar proof ap- 
peared in jTU] (see also [5]), and the proof provided here is based on similar ideas. 
We present the proof here because we require local existence for C 1 '" initial data 
(the proof in [10] used C 2,a ), and to ensure that bounds and existence time therein 
are independent of N. 

Without loss, we assume to — 0, and Uq, . . . , Uq to be !Fq measurable. Let U be 
a large constant and T a small time, both of which will be specified later. 

Define U — U(T,U) be the set of all time continuous ^-(-adapted C 1,a valued 
divergence free and spatially periodic processes u such that 



JV t~i te[o,r] 

hold almost surely. Also, we define M. = Ai(T) to be the set of all time continuous 
^-adapted C 1 ' 01 valued spatially periodic processes fi such that 



dXf u = u t {Xl' u ) dt + y/2^dWl 
with initial data X^ u {y) = y. As before, define F/ ,u = (A^'")" 1 , and define 

\r(y)=xr(y)-y 



Clearly a fixed point of W will produce a solution to the system l|3.ip - l|3.2p . Thus 
the proof will be complete if we show that for an appropriate choice of T and U, 
W maps Li into itself, and is a contraction with respect to the weaker norm 






»=i 



\\u\\ u = sup \\u t \\ a 
te[o,T] 

We first show W maps U into itself, using the two lemmas below. 
Lemma A.l. There exists c = c(a) such that 



||W(«)|| lia <c ma^l+IIV/Z'iJ + " ^Ell M o 



i=l 



a.s. 



Proof. First recall P vanishes on gradients. Thus 
(A.l) P[(V*Y)«] = -P[(V*u)y]. 



Now 



d t P [(V*Y)v\ = P [(V*Y)d t v + (V*<9;1>] 
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= P [{V*Y)diV - (V*v)diY] 

where we used (|A.1[) for the second term. Note that the right hand side involves only 
first order derivatives. Since P is a standard Calderon-Zygmund singular integral 
operator, which is bounded on Holder spaces, we obtain the estimate 

l|P[(v*y)«]lli,a<c||v*y|| a iH| liQ 

for some constant c = c(a). 

Applying this estimate to W, we have 

N 



(A.2) 



l|W(u)t|| lia < 



c 
N 



J2 \r Y t' 



a.s. 



from which the Lemma follows. 

Lemma A.2. There exists T = T(U, a) such that X i ' u ,n i ' u S M(T). 



□ 



We note that the diffusion coefficient is spatially constant, and thus we get the 
desired (almost sure) control on VA. Since V • u = 0, det(VAT l,tl ) = f , giving the 
desired control on V fj,. The details are standard, and we do not provide the them 
here (see for instance [TUI[T3"] ). 

Now choosing U = c(|) 2+Q -^ X)II u oIIi,oj an d then choosing T as in Lemma lA~2l 
Lemma \A . 1 1 shows that W maps U(U,T) into itself. Note that each summand on 
the right of (|A.2[) is bounded by U, which will prove the bound (|3.3[) . Further, 
given a uniform (in i) bound on ||uq||i, q , our choice of U can be made independent 
of N. 

It remains to show that W is a contraction. By definition of W we have 

iv 



W(u) t -W(v) t 



i=l 
N 



V*Y t hV 



u oY t 



iEp[(vT)(«5o^ 

i=l 

i=i 

1 n 

-^p[(v*rr)(<o^ 

i=i 

i N 

-^E p [ v *K°^)(n i,u -^ 



u oY t 



where we used the identity (|A. 1|) to obtain the last equality. Now we recall that 
jjL % ' u ,H % ' v £ M, and take C a norms. This gives 



(A.S) 



\W(u) t -W(v)t\\ a <^J2W u oh a 



N 



Y 



y: 



i=l 

Now from the definition of Y l ' u and Y % ' v we have 

ft 



y; 



y: 







[u s (y:> u )-v s (y:> v )] ds. 
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Taking C a norms, and applying GronwalPs inequality, and absorbing the exponen- 
tial in time factor into the constant c gives 

ft 

— vJ\„ ds a.s. 



JV 









a Jo 



Returning to (|A.3[) we have 

\\W(u) t - W(v) t \\ a < - sup K - «.|| i NIL a 



8<t 



Choosing t small enough one can ensure W is a contraction mapping. A standard 
iteration now shows the existence of a fixed point of W, concluding the proof of 
Proposition 13.11 



Appendix B. Logarithmic L°° bound on singular integral operators. 



In this appendix we provide a proof of Lemma 13.41 We restate it here for the 
readers convenience. 

Lemma (|3 .4(1 . If u is divergence free and periodic in Mr, then for any a € (0, 1), 
there exists a constant c = c(a, d) such that 



v«|| L oc ^ c H^H^o. 



i + hr 



where lj = V x u. 



Proof. Let K be a standard Calderon-Zygmund periodic kernel, and we define the 
operator T by 

Tf = K*f 

then we will prove that 

l/L 



(B.l) 



\\Tfh~ <c||/|| L . l + ln+ 



\f\\ 



This immediately implies Lemma 13.41 because we know Vu = — V(A _1 )V x u>, and 
— V(A _1 )Vx is a Calderon-Zygmund type singular integral operator. 

Now we prove (|B.1[) . We assume for convenience that all functions are periodic 
on the cube [0, l] d . We also recall that the kernel K satisfies the properties 

(1) K(y) < c\y\- d when \y\ < |. 

(2) Jj y|=p K(y) da{y) = for any r e (0, ±). 

Pick any e € (0, §). Then 



Tf(x) 



[0,1] 



K{y)f{x-y) dy 



(B.2) 



K (y)f(x - y)dy - 



K (y)f(x~y) d y 



\y\>£ 



Using property [T] about K, we bound the second integral by 



K{y)f{x-y) dy 



\v\>e 



<c||/|| 



<c ||/|| £ - In (i) 



„d-l 
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Using property [2] about K, we bound the first integral in (|B.2j) by 



f K(y)f(x-y) 



dy 



K(y)(f(x-y)-f(x))dy 



\y\<z 



<c\f\ a l \y\ a - d dy 

J\v\<e 

= c\f\ a s a 

Combining estimates we have 

\\Tf\\^^c[e a \f\ a + \\f\\ L ^ln^)} 

Choosing 

1 f\\fh~\ 1/a ' 



e = mm ■ 



l/l c 



finishes the proof. 



□ 
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